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ABSTRACT: 

in 

C\| ' Probability distributions of human displacements has been fit with ex- 

' ' c-| " ponentially truncated Levy flights or fat tailed Pareto inverse power law 

probability distributions. Thus, people usually stay within a given location 

O 

c/2 ■ (for example, the city of residence), but with a non- vanishing frequency they 

visit nearby or far locations too. Herein, we show that an important empiri- 
cal distribution of human displacements (range: from 1 to 1000 km) can be 

Oh' 

well fit by three consecutive Pareto distributions with simple integer expo- 
nents equal to 1, 2 and (<j) 3. These three exponents correspond to three 
displacement range zones of about 1 km < Ar < 10 km, 10 km < Ar < 

oo ■ 

300 km and 300 km < Ar < 1000 km, respectively. These three zones can 
be geographically and physically well determined as displacements within a 
city, visits to nearby cities that may occur within just one-day trips, and 



visit to far locations that may require multi-days trips. The incremental 
integer values of the three exponents can be easily explained with a three- 
scale mobility cost/benefit model for human displacements based on simple 
geometrical constrains. Essentially, people would divide the space into three 
major regions (close, medium and far distances) and would assume that the 
travel benefits are randomly /uniformly distributed mostly only within spe- 



cific urban-like areas. The three displacement distribution zones appear to 
be characterized by an integer (1, 2 or ^3) inverse power exponent because 
of the specific number (1, 2 or <;3) of cost mechanisms (each of which is 
proportional to the displacement length. The distributions in the first two 
zones would be associated to Pareto distributions with exponent f3 — 1 and 
(3 = 2 because of simple geometrical statistical considerations due to the 
a-priory assumptions that most benefits are searched in the urban area of 
the city of residence or in specific nearby cities. We also show, by using 
independent records of human mobility, that the proposed model predicts 
the statistical properties of human mobility below 1 km ranges, where peo- 
ple just walk. In the latter case, the threshold between zone 1 and zone 2 
may be around 100-200 meters and, perhaps, may have been evolutionary 
determined by the natural human high resolution visual range which char- 
acterize an area of interest where the benefits are assumed to be randomly 
and uniformly distributed. This rich and suggestive interpretation of human 
mobility may characterize other complex random walk phenomena that may 
also be described by a N-piece fit Pareto distributions with increasing inte- 
ger exponents. This study also suggests that distribution functions used to 
fit experimental probability distributions must be carefully chosen for not 
improperly obscure the physics underlying a physical phenomenon. 
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Human motion is a particular type of recurrent diffusion pro- 
cess where agents move from a location to another and return 
home after a given complex trajectory made of a certain number 
of steps of given lengths. For a large set of agents this motion can 
be considered equivalent to a complex random walk described by a 
peculiar statistics. Understanding the statistical nature of human 
motion is the topic of the present work. In the literature, human 
displacement distributions are traditionally fit using an exponen- 
tially truncated Pareto or inverse power law distribution based on a 
single inverse power-law exponent. However, the physical meaning 
of this specific distribution function, as well as of the measured in- 
verse power law exponent (for example, j3 = 1.75) have been never 
explained. Herein we show that the same human displacement 
distribution can be better fit with three inverse power law dis- 
tributions described by simple integer exponents (/3 = 1, 2, 3 or 
larger). Each exponent characterizes a specific displacement length 
zone whose ranges are about 1 km < Ar < 10 km, 10 km < Ar < 
300 km and 300 km < Ar < 1000, respectively. This alternative fit 
methodology for interpreting human displacement distributions is 
crucial because it suggests a clear physical/geometrical generating 
mechanism. Essentially, people would divide the space into three 
regions (close, medium and far distances), and they would assume 
that in each zone region the travel benefits are randomly /uniformly 
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distributed mostly within the urban area of the cities and would 
associate to each zone region an increasing integer number of cost 
mechanisms required to cover the given distance. We suggest that 
the probability associated to a displacement length is a monotonic 
function of its cost: less expensive (or shorter) displacements take 
place more likely than more expensive (or longer) displacements. 
Human mobility may be conditioned by multiple independent cost 
functions that work alone or may be statistically combined together 
according to the displacement length. For example, typical cost 
functions taken into account by people are the cost of the fuel, the 
time duration for the displacement, lodging costs that are required 
for very long trips, and others. The crucial fact is that on average 
these cost functions are likely directly proportional to the length 
of the displacement itself. For example, on average the fuel cost 
for covering a distance L is half the cost for covering a distance 2L: 
the same is true for the time interval needed to cover a give dis- 
tance. Thus, we can expect that the simplest strategy adopted by 
humans in optimizing their movements is that when fuel cost alone 
are taken into account by a person, displacements of length L are 
twice more probable than displacements of length 2L because of the 
homogeneous benefit distribution assumption within each urban 
zone. Thus, by assuming that people optimize their movement by 
uniformly distributing their travel resources (in time, energy and 
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money) within given range zones, we demonstrate by simple geo- 
metrical considerations based on area ratios that each displacement 
cost function would yield to a displacement probability distribution 
described by a basic inverse power law distribution P(Ar) oc 1/Ar. 
When nearby cities are visited a similar reasoning would yield to 
inverse power law with /3 — 2 because at leat 2 cost basic mecha- 
nisms would be activated and because the searched benefit would 
still be mostly searched within an urban zone. If 3 or more dis- 
placement cost functions condition the decision of an agent as for 
visits to far locations, then the combined probability distribution 
is the product of more basic probability distributions and would 
be characterized by inverse power law distributions with integer 
exponent (5 = 3 or more. We also show that a similar model pre- 
dicts the statistical properties of human displacements below 1 km 
ranges by people who just walk, and whose walking decision could 
be determined by time and physical energy considerations. In the 
case of just walking people the threshold between zone 1 and zone 
2 may be a few hundred meters and, perhaps, it may have been 
evolutionary determined by the natural human high resolution vi- 
sual range. The peculiarity of the three-scale pattern shown by the 
human displacement distribution suggests that on average humans 
take into account one, two, three (or more) alternative displace- 
ment cost mechanisms according to the length of the displacement 
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itself and assume that the benefits are uniformly /randomly located 
only within specific areas of interests, which characterize zone 1. 
This property may be more general and may characterize other 
complex random walk phenomena that are commonly interpreted 
in the scientific literature with simple exponentially truncated in- 
verse power law distributions based on a single exponent that may 
obscure, instead of clarifying, the physical mechanism behind a 
physical phenomenon. 

1 Introduction 

Complex random walks (from Brownian fractional diffusion to Levy walk 
and flights) are interesting physical phenomena that have found several ap- 
plications in physics |T] including human gait dynamics [2] . One of the most 
intriguing complex motion is the human one. Understanding human mobility 
laws is important for urban planning [3], traffic forecasting jUE], and for un- 
derstanding the spread of biological [HI [71 [8] and mobile viruses [9] . Human 
mobility also determines the formation of several social networks. In fact, 
in several cases the formation of human networks necessities the physical 
displacements of individuals. The social environment where humans move 
and make encounters, forming links and consequently networks, is mostly 
made of the residence city where people mostly live and work, and by a 
set of nearby or far locations and cities that can be occasionally visited. 
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Thus, understanding human mobility is also important for determining how 
specific social networks form and grow. Small-world social networks [10] as 
well as more complex networks such as free scale-social networks and others 
[TT] may be partially determined and explained by the specific patterns that 
characterize human mobility. 

Herein, we propose a specific model for understanding human mobility 
that takes into account the fact that a careful analysis of empirical data re- 
ferring to human displacements reveals that this phenomenon may be char- 
acterized by a multiple-scale Levy flight cost mechanism model constrained 
by simple spatial and temporal considerations that should constrain human 
mobility. 

The purpose of our study is to propose a possible physical explanation of 
human mobility interpreted in its overall statistical properties as phenomeno- 
logically manifested in probability distributions obtained from data collected 
from a very large sample of people. Thus, our study is not aimed to deter- 
mine the dynamics of the mobility of each person which, of course, will be 
very different from person to person and from location to location. How- 
ever, when the dynamics of a very large sample of persons in a multitude of 
locations is averaged in an overall probability distribution of displacements, 
the individual properties of each trajectory are smoother out and the emerg- 
ing statistical pattern would manifest a common underling characteristic of 
human mobility, which is the topic addressed in the present study. 

Empirical distributions referring to human and animal displacements have 
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been often approximated with various random walk or diffusion models that 
typically obey Pareto inverse power law or Levy flight distributions [121 H31 
dU [T51 [EE1 HZl [IB] . Inverse power law distributions are characterized by fat 
probability distribution tails. These tails indicate that while most human 
displacements are relatively short, the frequency of long displacements is not 
negligible either. Trips to nearby or far locations are rare, but not statistically 
impossible! If on average six links are sufficient to connect everybody to 
everybody else |T9], one explanation may be because of the peculiar inverse 
power law distribution that characterizes human mobility. In fact, the fat 
tail of these probability distributions implies that people can meet and form 
links also with individuals living quite far from their residence location. 

Pareto or inverse power law distributions are characterized by a specific 
power exponent j3 such as: 



Numerous physical, biological and social systems are characterized by Pareto 
distributions ^20\ . However, there are a few fundamental methodological and 
physical problems that need to be addressed for properly understanding the 
physics underlying a given complex phenomenon described by a probability 
distribution. 

1) Eq. ([I]) is a mathematical abstraction. No real system can be described 
by such a function for any real value of the variable x. In real systems the 
tail of the distribution is truncated because all natural systems are physically 




(1) 
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limited. 

2) Real complex systems are often not described by just one scaling expo- 
nent. Real complex systems are often multi-scales, that is, they are charac- 
terized by probability distributions that can be divided into different range 
zones. Each distribution zone may be characterized by a set of physical 
properties that determine the specific geometry of its probability distribu- 
tion within such a zone. Examples of bi- or multi-mode distributions are 
common [2Tj . 

3) Finally, the specific value of the measured exponents of these inverse 
power law distributions does have a physical origin that has to be properly 
investigated and explained. 

Probability distributions can also be fit with alternative functions gen- 
erated by very different physical background mechanisms, as the super- 
statistics of Beck and Cohen [22J has shown. In general, it may be even 
difficult to correctly determine which distribution function should be used to 
properly fit and interpret a given complex probability distribution. 

For example, some authors [23] have adopted exponentially truncated 
inverse power law function. Gonzalez et al. [12] and Edwards et al. [15] 
explicitly fit their human displacement distributions with an exponentially 
truncated inverse power law distribution because Eq. ([1]) would not fit the 
data. The same is done by Rhee et al. [H] where their displacement distri- 
butions obtained from people who mostly walk are fit with truncated Pareto 
distributions. The choice of these authors, however, is not fully backed by a 
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clear physical model and was proposed without a single explanation. 

Indeed, an exponential decay of the tail of a displacement distribution 
may be due to poor statistics or by geographical constrains that would not 
allow animals or humans to move for a distance longer than a certain maxi- 
mum. For example, it is reasonable to assume that the Levy flight patterns 
observed in wandering albatrosses may terminate with an exponential cut- 
off because these flights may be as long as several thousand kilometers and 
this geographical scale becomes compatible with the size of the Earth [13]. 
However, the issue may be more complicated and specific geometries of the 
displacement distribution may also indicate specific physical characteristics 
of human/animal movements that should be properly investigated. 

The purpose of this work is to study whether human displacement dis- 
tributions present more complex patterns that cannot be properly described 
by an exponentially truncated inverse power law distribution. In particular, 
we will show evidences that human displacement data can be described by a 
multi-scale inverse power law distribution that would eventually approach an 
exponential cut off. Then, we will propose that this multi-scale inverse power 
low distribution can be explained using a multi-scale statistical displacement 
cost mechanism grounded on commonsense physical and geometrical consid- 
erations that should characterize human mobility. 
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2 The traditional fit methodology: range from 
1 km to 1000 km 

One of the first large scale human motion studies was attempted by Brock- 
mann et al. [13]. These authors based their analysis of human movements 
on the trajectories of 464,670 dollar bills obtained from the bill-tracking 
system. From those reports these authors calculated the geographical dis- 
placements Ar = \xa — %b\ between a first (xa) and secondary (xb) report 
location of a bank note and the elapsed time AT between successive reports. 
They concluded that the distribution of these displacements obeys an inverse 
power law P(Ar) oc Ar _/3 with exponent /3 = 1.59. However, Brockmann 
et al.'s data do not allow to exactly track human movements, but, more 
properly, bank note movements that pass from hand to hand. The observed 
distribution does not reflect the motion of individual users but some kind of 
convolution between population-based heterogeneities and individual human 
trajectories. 

More recently, Gonzalez et al. [12J based their large scale study on trajec- 
tories of 100,000 anonymized mobile phone users. Each time a user initiated 
or received a call or a text message, the location of the mobile phone tower 
routeing the communication was recorded, allowing to reconstruct both user's 
time and space resolved trajectory within a reasonable resolution. In this way, 
these authors could measure the distance between user's positions at consec- 
utive calls. They collected a very large data set: 16,264,308 displacements 
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over a six-month period. This data set was labeled as 'Dl' in Gonzalez et al.. 
Because of its large size, this data set is probably the best available record 
for attempting to understand whether some simple statistical law governs 
human mobility. Herein we study the displacement distribution proposed by 
Gonzalez et al. that is depicted in Figure 1. 

The area covered by each mobile phone tower has a varying extension. 
Inside urban areas the distances between adjacent towers is lower than 1 
km. In the rural regions the mobile phone towers become sparser and their 
relative distance is up to a few kilometers. Thus, the distribution depicted in 
Figure 1 represents human mobility with the minimum resolution of about 
1 kilometer. The distribution range falls within 1000 km. This range is well 
below the geographical size of Europe where the data were collected from 
a cellphone provider. Thus, the entire range of the distribution for 1 km 
< Ar < 1000 km, or at least the range from 1 km to 300 km, likely reflects 
a real property of human mobility and should be explained by a physical 
model. 

Gonzalez et al. [12\ fit in their distribution (depicted in their figure 1C) 
with an exponentially truncated inverse power law (Levy flight) function of 
the type: 

P(Ar) oc (Ar + Ar ) _/3 exp(-Ar/K). (2) 

The fit gives: (3 = 1.75 ± 0.15, Atq = 1.5 km and cutoff values K rj 400 km. 
Analysis of the data also show that individual users return to a few highly 
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frequented locations with a daily frequency, such as home and work office, and 
there exists a weekly recurrence cycle as everybody would expect that humans 
do. This recurrent movement is quite different from, for example, bank notes 
that, instead, diffuse [13]. Thus, Gonzalezs record is more appropriate for 
studying human mobility. 

However, no physical explanation has been proposed by Gonzalez et al. 
to justify Eq. |2j It appears that the function fl2]) has been chosen to fit the 
data just because this is how this kind of distributions is traditionally fit in 
the literature, and because this function can approximately fit the data with 
appropriate regression parameters. No physical explanation of the found 
exponent = 1.75 has been proposed by Gonzalez et al. Thus, no real 
physical model has been proposed to explain the distribution depicted in 
Figure 1. 

3 An alternative fit methodology 

The present work focuses on an alternative interpretation of the distribution 
depicted in Figure 1. In fact, Figure 2 shows that the distribution of human 
displacement can be divided into three different zone ranges that can be fit 
with three inverse power laws of the type of Eq. ([I]): pp{r) oc l/Ar' 3 . We 
get: Zone 1, fa 1 for 1 < Ar < 10 km; Zone 2, fa 2 for 10 km < Ar < 300 
km; Zone 3, ~ 3 for 300 km < Ar < 1000 km. The distribution is assumed 
to be statistically meaningful within all represented ranges from about 1 km 
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to 1000 km, as claimed by Gonzalez et al. [12] . 

From a purely statistical perspective fitting the distribution with multiple 
inverse power law functions as done in Figure 2 appears to be statistically 
more accurate than fitting it as done in Figure 1 . This would be true at least 
within the range from 1 km to 300 km. In fact, the average \ 2 between the 
logarithm of the data and the fitting function model predictions gives: in the 
interval [1 km: 10 km], p = 0.0027 with Eq. © while p = 0.0020 with 
Eq. (J!) and = 1; in the interval [10 km: 300 km], ^ = 0.0039 with Eq. 
$2]) while = 0.0021 with Eq. (□) and = 2; in the interval [300 km: 1000 
km], p = 0.020 with Eq. (0) while x 5 = 0.026 with Eq. (pQ) and = 3. 

Thus, for displacements larger than 300 km, Eq. ([2]) performs slightly 
better than Eq. (JTJ), but this data are statistically poorer than those cov- 
ering the range from 1 km to 300 km, in particular at the very end of the 
distribution. Consequently, it may be possible the data may follow a Pareto 
law (JTJ) with = 3, as well as they may present an exponential cut off or, 
alternatively, they may be described with an inverse power law with a larger 
integer exponent. However, for displacements from 1 km to 300 km, where 
the data are more numerous and the distribution is statistically more robust, 
the distribution is clearly better fit with two consecutive inverse power law 
distributions with = 1 and = 2, which refer to two different displacement 
ranges. 

At a first sight the difference between the two alternative ways (depicted 
in Figure 1 and 2, respectively) of fitting the distribution of human displace- 
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ments appears to be so small that most people would consider the entire 
exercise to be irrelevant. However, a vigilant reader would recognize that 
the two alternative fitting methodologies are physically quite different. In 
the following section we show that the latter method (depicted in Figure 2) 
suggests an evident physical model of human mobility, while the former one 
(depicted in Figure 1) does not have any evident physical meaning. Thus, 
our exercise is not meaningless. It is a step forward a physical understanding 
of the complexity of the Levy-walk nature of human mobility based on a 
careful analysis of the data, which a simple exponentially truncated inverse 
power law distribution function with a single exponent might have obscured. 

4 A simple multi-scale diffusion model of hu- 
man displacements and its statistical/physical 
interpretation 

The interpretation of the above finding takes into account that geometrical 
constrains and travel costs both in terms of time and money determine the 
frequency of any kind of human displacement [21]. This behavior appears 
to be a general economical characteristic. In fact, in Krugman's economic 
geography [25], it is argued that people tend to cluster in cities (instead 
of distributing themselves homogeneously on the planet) partly because it 
saves on the transportation costs both in terms of time and money. The 
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clustering of the benefits in restricted areas produces geometrical constrains 
that force the distributions of displacements to assume an inverse power 
law distributions characterized by different integer exponents according to 
whether an agent is moving within an area or between area regions. 

I fact, according to the results depicted in Figure 2, the distribution of 
human displacements can be approximately separated into three consecu- 
tive displacement zones (relative to short, medium and long displacements). 
Each zone can be fit with three inverse power law functions with exponents 
approximately equal to 1, 2 and (~)3, respectively. This result suggests that 
human mobility may have a simple statistical explanation that can be based 
on simple displacement cost functions [24] and geographical topology. 

The almost sudden statistical transition observed from a displacement 
zone to the following one suggests that the number of cost functions involved 
in the process increases with the number that characterizes the displacement 
zone (1, 2 or 3). Thus, we propose a statistical model based on the empirical 
finding itself shown in figure 2 that assumes that the displacement probability 
distribution of human mobility is given by 



where the functions l//j(Ar) are the probability functions associated to in- 
dependent cost mechanisms that a given displacement Ar requires. Note 
that taking into account more than one cost mechanism should introduce 



iV(Ar) 



1 




(3) 
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a multiplicative term in the probability distribution because in statistics if 
A and B are two independent events the probability of their union is equal 
to the product of the two probabilities: P(AU B) = P(A)P(B). In our 
case it is evident that if two or more independent cost mechanisms are taken 
into account by an agent, the overall probability would be a weighted func- 
tion of all single probabilities associated to each cost function: that is, the 
overall probability is given by the product of the single probability functions 
associated to each cost mechanism, as well known in statistics. 

The peculiarity of the model is that the number of cost mechanisms in- 
volved in a displacement increases with the displacement zone number. Thus, 
we have: iV(Ar) = 1 for Ar < Ar x ; N(Ar) = 2 for Ar x < Ar < Ar 2 ; 
iV(Ar) = 3 for Ar 2 < Ar < Ar 3 , and potentially up to N zones. The 
displacements Arj with i = 1, 2, 3, . . . , N are the displacement borders that 
separate the zones. The model can be easily extended to N zones. 

Eq. [3] also implies that 

JV(Ar) 

P(Ar) J] /*( Ar ) = C N(Ar h (4) 

1=1 

where Cj with j = 1,2,3,... are specific constants that characterize each 
displacement zone. The constants implicit in Eq. H] would imply that a 
typical individual first divides the entire space around him in three major 
zones (close, medium and far distances). In fact, classifying things into three 
groups by size (small, medium, and large) is a quite natural methodology, as 
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also Aristotle would teach us. 

Eq. H] may also indicate that from a statistical perspective people opti- 
mize their movement by homogeneously distributing their resources in the 
geographic space available. This also would imply a travel length frequency 
inversely proportional to the overall (money and/or time) travel cost, which 
would yield to Eq. 

Costs that people take into account are travel time, fuel and lodging ex- 
penses, which are constrained by money budgets [26] . The crucial observation 
is that travel fuel cost and time duration are approximately proportional to 
the displacement length itself. Moreover, travel fuel cost and time duration 
are physically independent and different cost functions. Thus, for example, 
we can assume that if a person takes into account only the cost of the fuel 
or only the time duration for moving from a location A to a location B, a 
displacement of a given length Ar will be twice as more probable than a 
displacement of a given length 2Ar. This property, together with the as- 
sumption that the benefits are randomly and uniformly distributed within 
an urban area, should yield, for each cost mechanism, to Pareto distributions 
with exponent (3 = 1. 

The proposed hypothesis that the probability distribution associated to 
a single cost function is given by a Pareto distributions of the displacement 
length with exponent = 1 should be understood as a direct consequence 
of the observed Levy nature of the distribution of human displacements, as 
implicit in the probability distribution depicted in Figure 2. In fact, it is the 
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Levy nature of the distribution of human displacements itself that appears to 
exclude other possible distribution functions of the displacement such as, for 
example, an uniform distribution or an exponential distribution or a Poisson 
distribution or others, which would not be characterized by inverse power 
law Levy's tails. 

A simple statistical explanation of why within an urban area the proba- 
bility associated to a cost mechanism determining a given displacement may 
be well described by Pareto distributions of the displacement length with 
exponent (3 = 1 can be easily argued on geometrical considerations alone 
that would constrain the statistics by favoring the emergence of such dis- 
tributions. In fact, if the benefits are assumed to be uniformly/randomly 
distributed within an area of interest such as within a city, one of the reasons 
why a person may decide to move by a certain distance Ar = R is because 
s/he cannot find the desired benefit within a shorter distance Ar < R. Thus, 
two independent events must simultaneously occur: event A) the desired ben- 
efit is not found within the distance Ar < R: event B) the desired benefit 
is found at the distance Ar = R. Moreover, because we are assuming that 
within an urban area the benefits are uniformly/randomly distributed on a 
surface, the probability goes with the covered area; the probability of not 
finding a desired benefit within the distance Ar < R must decrease with the 
area of a circle of radius R and should be proportional to R~ 2 , while the 
probability of finding a benefit at a distance Ar = R must increase with the 
area of a ring of radius R and should be proportional to R itself. Thus, the 
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combined probability between event A and event B is given by the following 
equation: 

P(A UB) = P(A)P(B) oc x R oc 1 (5) 

which is a Pareto distributions of the displacement length with exponent 
(3 = 1. In conclusion, we suggest that the statistics associated to the above 
geometry forces the displacement distribution of human movement, which 
should be directly related to the displacement cost, to adopt an inverse power 
law formulation with exponent (3 = 1. Because this is the simplest case, we 
can assume that when a single cost mechanism is taken into account, it 
naturally produces Pareto displacement histograms with exponent = 1 as 
constrained by geometrical properties. 

Zone 1 refers to displacements from about 1 km to about 10 km that 
mostly refers to movements inside a city. The city where a individual is lo- 
cated is likely felt as a close zone. Within this range Figure 2 shows that 
the probability distribution is approximately proportional to 1/Ar. In fact, 
within an urban zone, a person may be directly influenced by the geometri- 
cal statistical constrains expressed in Eq. |5]and s/he may take into account 
only one cost mechanism. For example, the decision of moving by a certain 
distance Ar may be conditioned by fuel cost considerations or by time dura- 
tion considerations, but not by both. Thus, according to this interpretation, 
within the range of a few kilometers, each person thinks about saving fuel or, 



20 



alternatively, saving time, but on average only one cost mechanism propor- 
tional to the displacement distance Ar is taken into account by each person 
for each displacement decision, and the distribution associated to this cost 
would be constrained by the geometrical property expressed by Eq. |5j So, for 
short displacements we have that their frequency is inversely proportional to 
the travel cost (time duration or fuel consumption) which is approximately 
proportional to the travel length: 



P(Ar) ex — ]— ~ (6) 
D{Ar) Ar 



where D(Ar) is the physical cost function associated with the travel distance 
Ar. Consequently, ~ 1. 

Zone 2 refers to displacements larger than 10 km and shorter than 300 
km. These displacements refer mostly to short visits to nearby cities and a 
typical individual would feel this region as a medium zone. In fact, a person 
can cover such distances up to a few hours with a car or train and round 
trips within this zone may occur just within the day. Within this range 
the displacement probability distribution is approximately proportional to 
1/Ar 2 , as depicted in Figure 2. This pattern may emerge because within 
this displacement range both fuel cost and time duration costs may be taken 
in serious consideration by a traveler, who would prefer to save both fuel 
and time. Note that for displacements between 10 km and 300 km both 
the displacement time duration and fuel costs are not negligible for most 
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people. Thus, for displacements between 10 km and 300 km, that is for trips 
to nearby cities and towns, the displacement probability is likely the product 
of two probability functions which are both proportional to 1/Ar. For these 
displacements we have 

P(Ar) K dW) t(a7) ~ (7) 

where T(Ar) is the temporal cost function associated with the displacement. 
Consequently, the resulting inverse power law exponent is /3 ~ 2. 

The emergence of the above statistical pattern with an inverse power law 
exponent with (3 = 2 may also be favored by a geometrical constrain similar 
to that expressed in Eq. |5j In fact, it may be that each individual is looking 
for a rare benefit not found in his/her city of residence. An individual may 
feel that the desired rare benefit is located just in one specific nearby cities 
and would a-priory disregard other possible nearby cities as well as most of 
the country area region among cities, which would a-priory be considered 
irrelevant. In this case, while the probability of not finding a desired benefit 
within the distance Ar < R would still decrease with the area of a circle 
of radius R and should be proportional to R~ 2 , the probability of finding a 
benefit at a distance Ar = R would be likely proportional only to the size 
of a specific nearby city itself and, therefore, it would be independent of the 
distance R between the two cities. In this specific situation the combined 
probability between event A and event B is given by the following equation: 
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P(A U B) = P{A)P{B) oc — x cityarea oc — , (8) 

R 2 R 2 

which is a Pareto distributions of the displacement length with exponent 
(3 = 2. Thus, for travels to nearby cities an agent may assume that the desired 
benefits are not uniformly distributed within the entire available territory 
area, but s/he appears to be guided by a pre-decision that limits the search 
inside a restricted region which is more directly related to the area of a 
specific nearby city than to the distance between the residence city and the 
visited city. 

Zone 3 refers to displacements larger than 300 km. These displacements 
refer mostly to trips to far cities which are likely felt by individuals as a 
far zone because they require multiple day trips. In the eventuality that 
the data were collected from a phone provider localized in one of the minor 
European country, let us say Hungary which has an extension of 200-500 km, 
these far displacements would mostly correspond to abroad trips. Within this 
range, Figure 2 shows that the displacement probability distribution may be 
characterized by an exponent (3 = 3 and may end with an higher exponent 
or an exponential cut off. In this situation, the observed pattern may also be 
due to poor statistics or to the fact that for displacements longer than 300 
km some geographical boundary is approached. 

However, this third zone may be approximately fit with an equation pro- 
portional to 1/Ar 3 suggesting that at least one additional cost mechanism 
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proportional to the displacement length may be activated for these long trips. 
For example, in addition to travel fuel and duration costs, a traveler would 
have lodging costs because s/he will not be able to return home the same 
day. These additional costs may also increase with the displacement length 
because a traveler would try to optimize his economical resources by staying 
longer in a farer location for better exploring it by taking into advantage the 
rare visit. Thus, we can expect that the displacement probability distribu- 
tion relative to this third displacement zone may decay as an inverse power 
law with exponent ft = 3 or a higher order, and may eventually end with an 
exponential cut off. So, we have: 



P(M oc ^ T H(Ar) ~ Ar 3 ' (9) 



where H(Ar) ~ Ar + 0(Ar 2 ) is the cost function associated to the additional 
costs that these long trips require. Consequently, the inverse power law 
exponent is ft <J 3. 

Figure 3 depicts a schematic representation of human motion with Levy 
flight displacements in the three different displacement zone ranges. The 
shadowed circles indicate the urban area of the cities where most displace- 
ments occur. Within the shadowed city area the benefits are likely uni- 
formly/randomly distributed and the displacement distribution obey to a 
Pareto law with exponent ft — 1 as predicted by Eq. |5j When people 
move to nearby cities the searched benefits are still assumed to be uni- 
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formly /randomly distributed mostly within the shadowed cities zone, and 
the displacement distribution obey to a Pareto law with exponent /3 = 2 as 
predicted by Eq. [SJ When far cities are visited, the displacement distribution 
obey to a Pareto law with exponent /3 ~ 3 because of additional activated 
cost mechanisms. 

5 Displacement distributions of walking peo- 
ple: range below 1 km 

An interesting exercise would be to determine whether the above model may 
somehow be extended to describe and/or predict the behavior of people who 
just walk below the 1 km region. The data we have analyzed in Sections 2 
and 3 do not have such a resolution and refer to displacements from 1 km to 
1000 km obtained by means of standard vehicles such buses, trains 

and airplanes. Adopting displacements larger than 1 km filters out possible 
influences due to the small-scale topology of the local geography. Moreover, 
the fact that standards equivalent transportation vehicles are used by all 
agents implies that the travel costs both in money and time are evaluated 
equally by all agents. Consequently, those costs are in first approximation 
only function of the displacement length and would not significantly depend 
on the physical characteristic of an individual person. In fact, in these cases, 
the costs in money and/or time for covering an equal displacement would be 
approximately the same for a man and a woman, for a young person and 

25 



for an elderly, for an athletic person and a sedentary one. These common 
properties would favor the emergence of clearer distribution exponents very 
close to integer numbers such as (3 = 1 and /3 = 2, which would depend only 
on the displacement range zone, as explained above. 

About the human mobility behavior below the 1 km region, where people 
more likely walk, we still expect that the probability distributions of their 
displacements is determined by inverse power laws because of the general 
theoretical rule discussed in the previous section. In these cases, however, 
the subjective physical characteristics of the walker, as well as the small- 
scale topology of the local geography may induce a large variation among 
the detected exponents. However, because the major cost functions would 
still depend on one or, at most, two major cost components such as time and 
physical energy required for each displacement, we would expect that the 
behavior of walkers below the 1 km region may theoretically mirror at a lower 
scale the statistical behavior observed at larger scales. Also the geometrical 
conditions expressed in Eqs. [5] and M would be mirrored at a lower scales 
because people may decide to randomly move only within specific limited 
area regions where Eq. [5] works and between these area regions where Eq. [8] 
works. In other words, we would expect that the adoption of vehicles such 
clS CclXS and trains only magnify the range of human displacement, but does 
change the general properties that humans adopt to decide the frequency of 
their displacements and how they organize the space around them. 

Thus, we do expect that also when people just walk, they still divide the 



26 



space into three region zones (close, medium, far distances) as depicted in a 
larger scale in Figure 2, and that the displacement distributions of walking 
people may be characterized by exponents j3 ~ 1 at the smallest scale end (for 
example Ar < 100-200 m) and (3 « 2 at larger scale (100-200 < Ar < 1000 
m). 

Note that the far zone (above about 1 km) will unlikely be covered by 
simply walking by most people in our days, and distribution of experimental 
data not appropriately designed to exclude automated vehicles (for example, 
cars and buses) would not clearly manifest the third zone. In fact, if auto- 
mated vehicles were used for the covering these larger distances, the statistics 
implemented in Figure 2 would be activated. Thus, typical distribution of 
mostly walking people would at most manifest two zones. 

The expected 100-200 meter threshold between displacement zone 1 and 
displacement zone 2, which are characterized by just one or two activated 
cost mechanisms, respectively, may be associated to the typical human high 
resolution visual range area, which limit a natural range around each person, 
and may have been determined by evolution. In fact, at smaller scales (less 
than 100-200 meters) a typical person would be able to see even small details 
such as small objects and animals. This visual range would constitute a 
natural close range and would be psychologically equivalent to the residence 
city space once that fast vehicles are used to cover these larger distances 
within just a few minutes. In this close range zone the decision of moving may 
be determined by only one cost mechanism and it would be treated as an area 
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where the benefits are randomly/uniformely distributed. This configuration 
activates the statistics implemented in Eq. [5j Thus, we should expect Pareto 
distributions with (3 ~ 1. On the contrary, for longer scales (larger than 100- 
200 meters) the displacement has to be more carefully planned because it 
would involve locations not in the immediate vicinity of the agent (like in 
the case of visiting nearby cities and towns), and we should expect that the 
displacement decision is made of two cost mechanisms and the geometrical 
configuration of displacements between areas of interests activates Eq. [HJ 
Consequently, we should expect Pareto distributions with /3 ~ 2. 

The above exponents, (3 = 1 and (3 = 2 respectively, should be expected 
in the cases in which the individual movement is not strongly influenced by 
specific environmental topology. In fact, when people are forced by environ- 
mental topology to cover larger distances more often and small distanced less 
often, we would expect approximate Pareto distributions with 1 < (3 < 2, 
where the threshold between zone 1 and zone 2 becomes less evident. 

The above predictions appear to be fully compatible with the recent re- 
sults of Rhee et al. [18] (see their figures 5 and 6), who studied the Levy- walk 
nature of human mobility of mostly walking people based on 226 daily GPS 
traces collected from 101 volunteers in five different outdoor sites such as two 
university campuses, in Manhattan (NY), in Disney World in Orlando (FL), 
and one state fair. 

Rhee et al. [18] used truncated Pareto law to fit his distributions which 
may present the same problem of occulting the real physics of the process. 
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Interestingly, among the five records studied by Rhee et al., the one that 
detected an inverse power law exponent closer to 1 ((3 = 1.16 within the range 
0.1-10 km), which would correspond to our finding of figure 2, was the one 
relative to Manhattan (NY) where the displacement means included vehicles 
such as cars and buses, as for the data used for obtaining the distribution 
depicted in Figures 1 and 2. Movements inside a city also more likely activates 
the statistics implemented in Eq. |5j Thus, we should have expected a Pareto 
distribution with 

We observe that the walking in the two university campuses (North Car- 
olina State University (3 = 1.29; Korea Advanced Institute of Science and 
Technology (3 = 1.66) could be strongly constrained by geographical topol- 
ogy because students are more likely forced to move to specific locations, 
and here we found Pareto distributions with 1 < < 2, where the threshold 
between zone 1 and zone 2 is not evident. 

On the contrary, the walking in the Disney World (/3 = 1.82) and the 
state fair in Raleigh (NC) (/3 = 1.74) is likely less constrained by specific 
geographical topological necessities, and these data may more accurately de- 
scribe general properties of free human mobility. However, it is an impression 
of the present author that the distributions depicted in figure 6d and 6e in 
Rhee et al. [18] referring to the latter two cases would be better fit with a 
two-scale Pareto law, and present a threshold from a low Pareto exponent 
closer to /3 — 1 to a high Pareto exponent closer to (3 = 2 around 100-200 
m. This claim is proven herein in Figure 4 that reproduces figure 6d in Rhee 
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et al., where the pattern is clearer because apparently obtained with more 
data. 

The finding depicted in Figure 4 appears to confirm our expectation that 
the behavior of simple walkers below the 1 km region approximately mirror 
at a lower scale the behavior observed at larger scales depicted in Figure 2. 
We see two displacement zones characterized by /3 = 1 and (3 = 2, respec- 
tively, where one or two cost mechanisms are likely taken into account by 
a walker according to whether the displacement is felt as close or medium. 
Moreover, we notice that Disney World park is organized e;e;lomer- 
ate of little town/parks that resemble in a miniaturized scale the typical 
geographical topology of a set of close urban areas. Thus, it may be not sur- 
prising that the displacement distribution depicted in Figure 4 looks similar 
to that depicted in Figure 2 by reproducing at a smaller scale the geometri- 
cal/statistical properties of Zone 1 and Zone 2 that activate the conditions 
expressed in Eqs. [5] and [HJ which yield to Pareto distributions with (3 = 1 
and (3 = 2, respectively. 

6 Conclusion 

In this work we have studied the distribution depicted in Figure 1 that is 
made of a sample of several million of human displacement data [12] . These 
authors have traditionally fit their distribution with a typical exponentially 
truncated inverse power law with one exponent without being able to explain 
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its physical origin nor the physical meaning of the measured exponent. 

However, it is legitimate to look for alternative ways of fitting a complex 
probability distribution [22] if this exercise yields to a physical explanation of 
the phenomenon under study. We have found that the same distribution may 
be fit with a multi-scale inverse power law distribution with simple integer 
exponents 1, 2 and 3 within three zones that depend on the length of the 
displacement: from 1 to 10 km, from 10 to 300 km and from 300 to 1000 km, 
respectively. 

The multiple integer exponent distribution may be generated by the fact 
that people optimize their movements by first dividing the entire available 
space into three zones corresponding to close, medium and far distances. 
Then they assume that on average the benefits of travel are distributed on 
average randomly and uniformly within specific area of interests such as cities 
at all travel distances. Finally, people would take into consideration a incre- 
mental number of cost mechanisms according to each zone, which are always 
approximately proportional to the displacement length such as travel fuel 
cost, travel time duration and others. Thus, as we have proven, each sin- 
gle mechanism alone would yield to a Pareto probability distribution with 
exponent — 1 because this kind of distribution would fit a geometrical 
constrain related to the condition of benefits randomly and uniformly within 
specific area of interests as expressed in Eq. [5] that yields to Pareto distri- 
bution with (3 = 1. The combination of two or more cost mechanisms would 
yield to Pareto probability distribution with exponents j3 = 1, 2, 3 and so 
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on, according to the displacement zone that activates 1, 2, 3 (and so on) cost 
mechanisms. In particular, trips to nearby cities would be constrained by a 
geometrical condition expressed in Eq. |H]that yields to Pareto distribution 
with — 2. This geometrical constrain requires that when a person search a 
benefit in a nearby city s/he would still search it within a specific city/area 
like region, and s/he would a-priory disregard any other region at equivalent 
distance in other directions. 

The above conclusion appears banal and predictable. Everybody knows 
that when people travel they assume that specific benefits may be ran- 
domly/uniformly distributed in areas of interests such as within cities, and 
take into account things such as fuel cost, time duration etc. It is perfectly 
reasonable that a-priory, on average, individuals divide the space in three 
major zones (close, medium, far distances). Thus, because of simple geomet- 
rical constrains we should have expected that human displacements could 
be described by simple inverse power law distributions with simple integer 
exponents. However, the importance of our work is to have proved this fact 
for the first time by using an actual human mobility distribution. 

In Section 5 we have show that the same model is capable to predict and 
explain also the complexity of the Levy-walk nature of human mobility for 
mostly walking people, as Figure 4 shows. In fact, the adoption of vehicles 
simply magnifies the geographical space taken into account by an individual, 
but the physical mechanisms remain the same. Walking people, too, would 
first divide the space into areas of interests and three potential regions (close, 



32 



medium and far distances), would assume that in each area the benefits are 
homogeneously distributed and would associate to each region an increasing 
integer number of cost mechanisms. This yields to multi-scale Pareto distri- 
butions with incremental exponents /3 = 1, 2 (and so on). However, the third 
region may not appear in current experimental distributions for the reasons 
explained in the text. 

The multi-scale mobility model described in Eqs. [31 IU [5] and |8] can be 
quite general, indeed, and can efficiently explain several real mobility and 
diffusion processes where the displacement is reasonably a function of geo- 
metrical constrains and the number of cost functions involved in the process. 
This cost function geometrically increases with the displacement length every 
time that a certain displacement zone is crossed. This yields to apparently 
truncated power-law distributions that, indeed, are made of different power 
law distributions with incremental exponents (in the human mobility case 
(3 = f, 2, 3, . . .) referring to different displacement length ranges, as suggested 
by the data depicted in Figures 2 and 4. 

References 

[1] N. Scafetta, Fractal and Diffusion Entropy Analysis of Time Series: 
Theory, concepts, applications and computer codes for studying frac- 
tal noises and Lvy walk signals. (VDM Verlag Dr. Mller 2010). 



33 



[2] N. Scafetta, D. Marchi and B. J. West, "Understanding the complexity 
of human gait dynamics," Chaos 19, 026108 (2009). 

[3] M. W. Horner, and M. E. S O'Kelly, "Embedding economies of scale 
concepts for hub networks design," J. Transp. Geogr. 9, 255 (2001). 

[4] R. Kitamura, S. Fujii and E. Pas, "Time-use data, analysis and model- 
ing: toward the next generation of transportation planning methodolo- 
gies," Transport Policy 4, 225 (1997).. 

[5] R. Kitamura, C. Chen, R. M. Pendyala, and R. Narayaran, "Micro- 
simulation of daily activity-travel patterns for travel demand forecast- 
ing," Transportation 27, 25 (2000). 

[6] V. Colizza, A. Barrat, M. Barthelemy, A. -J. Valleron, and A. Vespignani, 
"Modeling the worldwide spread of pandemic influenza: baseline case 
and containment interventions," PLoS Medicine 4, 95 (2007). 

[7] S. Eubank et al., "Controlling epidemics in realistic urban social net- 
works," Nature 429, 180 (2004). 

[8] L. Hufnagel, D. Brockmann, and T. Geisel, "Forecast and control of 
epidemics in a globalized world," Proc. Natl Acad. Sci. USA 101, 15124 
(2004). 

[9] J. Kleinberg, "The wireless epidemic," Nature 449, 287 (2007). 



34 



[10] S. Schnettler, "A structured overview of 50 years of small- world re- 
search," Social Networks 31, 165 (2009). 

[11] R. Toivonen, L. Kovanen, M. Kivela, J. Onnela, J. Saramaki, and K. 
Kaski, "A comparative study of social network models: Network evo- 
lution models and nodal attribute models," Social Networks 31, 240 
(2009). 

[12] M. C. Gonzalez, C. A. Hidalgo, and A.-L. Barabasi, "Understanding 
individual human mobility patterns," Nature 453, 779 (2008). 

[13] D. D. Brockmann, L. Hufnagel, and T. Geisel, "The scaling laws of 
human travel," Nature 439, 462 (2006). 

[14] G. M. Viswanathan, et al. , "Levy flight search patterns of wandering 
albatrosses," Nature 381, 413 (1996). 

[15] A. M. Edwards et al., "Revisiting Levy flight search patterns of wander- 
ing albatrosses, bumblebees and deer," Nature 449, 1044 (2007). 

[16] G. Ramos- Fernandez et al., "Levy walk patterns in the foraging move- 
ments of spider monkeys (Ateles geoffroyi)," Behav. Ecol. Sociobiol. 273, 
1743 (2004). 

[17] D. W. Sims et al., "Scaling laws of marine predator search behaviour," 
Nature 451, 1098 (2008). 



35 



[18] I. Rhee, M. Shin, S. Hong, K. Lee, S.J. Kim, and S. Chong, "On the 
Levy- Walk Nature of Human Mobility," IEEE-ACM Transactions on 
Networking 19:3, 630 (2011). 

[19] D. J. Watts, Six Degrees: The Science of a Connected Age. (Norton & 
Company 2003). 

[20] M. Schroeder, Fractals, chaos, power laws: minutes from an infinite 
paradise. (Freeman, New York 1991). 

[21] M. Latapy, C. Magnien, and N. Del Vecchio, "Basic notions for the 
analysis of large two-mode networks," Social Networks 30, 31 (2008). 

[22] C. Beck and E. G. D. Cohen, "Superstatistics," Physica A 322, 267 
(2003). 

[23] T. Fenner, M. Levene, and G. Loizou, "A model for collaboration net- 
works giving rise to a power-law distribution with an exponential cutoff," 
Social Networks 29, 70 (2007). 

[24] E. H. Lundevaller, "The effect of travel cost on frequencies of shopping 
and recreational trips in Sweden," J. of Transport Geography 17, 208 
(2009). 

[25] M. Fujita, P. Krugman, and A. J. Venables, The Spatial Economy: 
Cities, Regions, and International Trade, The MIT Press (September 
1, 2001). 



36 



[26] P. Mokhtarian, and C. Chen, "Ttb or not ttb, that is the question: 
a review and analysis of the empirical literature on travel time (and 
money) budgets," Transportation Research Part A 38, 643 (2004). 



37 



0.001 



0.0001 



1e-005 



1e-006 







1.58* 




Displacement distribution 

(x+1 .5) A (-1 .75)*exp(-x/400) ■ 


o 
















































i \ 



1 10 100 1000 

Ar (Km) 

Figure 1: Distribution of human displacements (from record Dl and Figure 
1C in Gonzalez et al. [12]) • The length Ar is the travel distance. The 
distribution is fit with Eq. (J2J). 
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Figure 2: Distribution of human displacements (from record Dl and Figure 
1C in Gonzalez et al. |12j). The length Ar is the travel distance. The dis- 
tribution is approximately reproduced within three consecutive ranges with 
incremental inverse power law functions proportional to 1/Ar^ 3 . We get: 1) 
« 1 for 1 km < Ar < 10 km; 2) « 2 for 10 km < Ar < 300 km; 3) « 3 
for 300 km < Ar < 1000 km. 
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Figure 3: Schematic representation of human motion with Levy flights in 
three different displacement zone ranges. The bottom left corner represents 
the residence location where the walker returns. The regions with a higher 
density of displacements represent cities or towns (zone 1) where the walker 
stays longer. The distance among the cities classifies zone 2 or zone 3 accord- 
ing to whether the journey is medium or long. The shadowed areas represent 
the urban area where the benefits are assumed to be randomly /uniformly 
distributed. Displacements within these restricted areas yield to Pareto dis- 
tribution with exponent (3 = 1 according to Eq. [5j Displacements between 
close areas of interest yield to Pareto distribution with exponent = 2 ac- 
cording to Eq. [HI 
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Figure 4: Reproduction of figure 6 in Rhee et al. [18] which refers to a 
distribution of displacements at Disney World (Orlando, FL) by a sample of 
traces obtained from 18 volunteers that mainly walked in the park during 
holidays. The figure shows in green the authors' original fit obtained with a 
truncated Pareto distribution that gives an exponent (3 = 1.86 for the entire 
range. The red and blue curves, instead, were added by me and highlight 
that the smallest displacement range up to 100-200 meter is compatible with 
a Pareto distribution with (3 = 1, while the longer scale is compatible with 
a Pareto distribution with (3 = 2. Note that the statistics for displacements 
above approximately 1000 m is visibly poorer. 
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